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PROLOGUE 


A one dimensional, coupled ecosystem/carbon cycle model is used to analyze the 
biogeochemical-physical interactions and carbon fluxes in the Bermuda Atlantic Time-series 
Study (BATS) site for the period of 1992-1998. The influence of mesoscale eddies (remote 
forcing) is taken into account via nutricline displacements and vertical advection derived from 
satellite altimetry. Physical forcing and biogeochemical boundary conditions are derived from 
the comprehensive BATS data sets. The observed parameters not used for model forcing and 
boundary conditions are used to verify model performance. The model results compare well with 
observations (most variables are within 8% of observed values). The new production predicted 
by the model, 0.47±0.14 mol N m' 2 yr' 1 , is within 6% of the 0.50+0.15 mol N m yr mean 
estimate for the Sargasso Sea reported by Siegel et al. [1999]. The mean (1992-1998) sea-air CO 2 
flux was -0.3 mol C m' 2 yr' 1 , indicating that the BATS region is a weak sink of atmospheric C0 2 . 
This estimate was obtained with the Wanninkhof [1992] gas transfer algorithm, which provides a 
ApC0 2 range within 2.5% of the observed value. However, the sea-air flux varies depending 
upon the gas transfer algorithm used. Four of the five algorithms used provide a reasonable 
(<10%) match with the observed ApC0 2 . The sea-air flux obtained with the use of these 
algorithms ranges from -0.32 to -0.50 mol C m" yr 1 , which is within the range (-0.22 to -0.83 
mol C m' 2 yr' 1 ) of previously reported values [ Bates et al., 1996; Bates et al., 1998d] . The model 
reveals prominent seasonal TC0 2 variability, primarily due to SST and total alkalinity (TA) 
changes. The seasonal variability of surface layer TA is quite large. For the period of model 
simulation (1992-1998), the TA ranged between a low of 2364 pmol kg' 1 and a high of 2404 
|i.mol kg' 1 . A significant portion of the TA variability is due to salinity changes via the effects of 
evaporation and precipitation [Bates et al., 1996]. The minimum TA (2364 pmol kg ) occurred 
on August 1998, during the peak of La Nina in the tropical Pacific. This minimum TA occurred 
during the period’s warmest SST (28.4°C), lowest surface salinity (36.23), highest pC0 2 (420 
p.atm), and lowest TC0 2 (2018 nmol kg' 1 ). The overall carbon balance (fluxes through the 
surface and bottom boundary at 350 meters) consists atmospheric C0 2 uptake of 0.3 mol C m 
yr' 1 , upward DIC bottom flux via mixing and advection of 1.1 mol C m 2 yr ', and carbon export 
of 1.4 mol C m’ 2 yr’ 1 via sedimentation. Upper ocean TC0 2 levels increased between 1992 and 
1996 at a rate of -1.2 p.mol kg' 1 yr' 1 , consistent with observations [Bates et al., 1996; Bates, 
2001]. However, this trend was reversed during 1997-1998 to -2.7 pmol kg 1 yr *. This large 
change in TC0 2 trend coincided with a significant TA drop (-9 pmol kg' 1 ) from 1997 to 1998, in 
response to hydrographic changes imposed by the El Nino-La Nina transition, which were 
manifested in the Sargasso Sea by the warmest SST and lowest surface salinity of the period 
(1992-1998). 
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1. Introduction 

Bermuda is the site of some of the longest-running 
ocean and atmosphere time-series investigations. The 
linchpin of oceanographic time-series in Bermuda is the 
Hydrostation S program, started by Henry Stommel 
(Woods Hole Oceanographic Institution) and 
co-workers in 1954. Hydrostation S is located -22 km 
southeast of Bermuda (32° 10'N, 64° 30’W). The 
comprehensive long-term physical and biogeochemical 
data originating from Hydrostation S are invaluable for 
interpreting the interannual patterns in the physical 
structure of the ocean [e.g., Schroeder et al., 1959; 
Schroeder and Stommel 1969; Pocklington, 1972; 
Wunsch, 1972], the seasonal cycle of nutrients, 
chlorophyll, and primary production in the Sargasso Sea 
[e.g., Menzel and Ryther , 1960; Menzel and Ryther 
,1961], and ocean biogeochemistry in general. 

In October 1988 the Bermuda Atlantic Time-series 
Study (BATS) commenced sampling the Sargasso Sea 
in an area 85 km southeast of Bermuda (31°45'N, 
64°10’W) as part of the U.S. Joint Global Ocean Flux 
Study (JGOFS). The scientific goal of the BATS 
program is to understand the causes of seasonal and 
interannual variability in ocean biogeochemistry, both 
at the BATS study site and as it may relate to 
biogeochemistry of the rest of the ocean [ Michaels and 
Knap , 1996]. 

The biogeochemistry near Bermuda is influenced 
by strong seasonal meridional gradients. There is weak 
mixing to the south of Bermuda which leads to a 
permanently stratified water column of oligotrophic 
characteristics throughout the year [Michaels and Knap , 
1996]. Nitrate is rarely present at the surface, the 
ecosystem is dominated by picoplankton and the 
microbial food web, and nitrogen-fixing cyanobacteria 
are common [ Carpenter , 1983; Olson et al., 1990a, b; 
Siegel et al., 1990]. In contrast, the deep winter mixed 
layers to the north of Bermuda result in blooms of 
larger phytoplankton, including diatoms and 
coccolithophores, and a complex transition to 
oligotrophy in the summer [ Siegel et al, 1990]. Both of 
these seasonal patterns can occur at the BATS site 
depending on the intensity of winter mixing [Michaels 
and Knap , 1996]. 

Many unanswered questions remain regarding the 
nutrient and carbon balances at BATS. Previous efforts 
to achieve closure on the carbon budget at BATS 
[Michaels et al., 1994b] revealed a large residual 
discrepancy (-8 mmol C m 2 d' 1 ). The uncertainties in 
the vertical fluxes were very small except for two 
processes, particle settling and horizontal advection. A 
first-order carbon balance required either that the 


shallow sediment traps under-collected sinking particles 
by as much as a factor of 5 or that up to 70% of the 
carbon dynamics are determined by advection, or some 
combination of both [Michaels et al., 1994b]. 

In this study, we combine existing long-term data 
sets with the use of a one-dimensional ecosystem model 
to investigate the physical and biogeochemical 
processes that govern the carbon and nutrient cycles at 
BATS. Considering the limitations of the 

one-dimensional dynamics, the influence ot mesoscale 
eddies (remote forcing) is taken into account via 
nutricline displacements and vertical advection derived 
from satellite altimetry. The outstanding question of 
how influential horizontal advection is to the seasonal 
fluxes of nutrients and carbon will remain unanswerable 
until three-dimensional field studies are conducted at 
BATS. However, our study shows that a large 

uncertainty remains in the determination of the air-sea 
C0 2 flux depending on the formulation used. One of 
the main objectives of this study is to evaluate the 
effects of the different gas exchange formulations on 
the carbon balance in the Sargasso Sea. In addition, the 
model provides an extension for the seasonal and 
interannual analysis beyond the available observations. 
This is more evident for the ApC0 2 record, for which 
only 2.5 years of observations are currently available 
(1994-1996). 

2. Data Sources 

The ecosystem model is physically forced by 
temperature, salinity, and mixed layer depth derived 
from BATS data. The data sets from BATS are 
discussed in several previous studies [Bates et al, 1996; 
Michaels et al., 1994a; Michaels and Knap, 1996]. The 
BATS biogeochemical data were used for model 
validation. The mixed layer and PAR time series were 
obtained from the Bermuda Bio-Optics Program 
(BBOP) data base. The mixed layer depth (D,„| d ) was 
calculated by finding the first depth where <7e(D m i d ) - 
ct 9 (0) = a T AT; a T is the coefficient of thermal 
expansion at sea surface conditions and AT is chosen to 
be 0.5°C [Siegel et al., 1995]. 

The 1994-1996 seawater and atmospheric pC0 2 data 
originate from the carbon dioxide measurements 
obtained from the R/V Weatherbird II in the Sargasso 
Sea [Bates et al., 1998a; Bates et al., 1998b; Bates et 
al., 1998c]. Most of these data sets were obtained 
interactively via the web-based data extractor 
maintained by BBSR (Bermuda Biological Station for 
Research, Inc.; www.bbsr.edu/cintoo/bats/bats.html). 
The eddy-induced vertical velocity time series was 
obtained from the divergence of the mesoscale eddy 
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velocity field derived geostrophically from the merged 
TOPEX/Poseidon-European Remote Sensing Satellite 
(ERS), 0.25°, 10-day sea-surface height anomaly 
(SSHA) grids [AVISO, 1997], The SSHA data were 
obtained via anonymous ftp from the AVISO data 
distribution server (ftp://ftp.cls.fr/pub/oceano/AVISO). 
The wind speed time series required for the gas transfer 
formulation was obtained from the Bermuda Airport 
weather station. The gridded winds to calculate the 
Ekman pumping vertical velocity, and cloud cover data 
to calculate surface PAR, were obtained from the 
NCEP/NCAR Reanalysis Project [Kalnay et al., 1996] 
daily products, extracted from the Lamont-Doherty 
Earth Observatory (LDEO) web-based data library 
(ingrid.ldeo.columbia.edu/SOURCES/). 

3. Model Description 

The model evolved from a four-component 
ecosystem model [McClain et ai , 1996] applied to the 
study of physical and biological processes at Ocean 
Weather Station P (OWS P). There have been three 
upgrading stages of development for the model since 
the McClain et al. [1996] version. The first 
enhancement was conducted for the equatorial Pacific 
Warm Pool, which included a spectral formulation for 
the downwelling irradiance, nitrification of NUf to 
N0 3 , and allowed physical forcing from offline 
coupling to a three-dimensional ocean circulation model 
simulation [McClain et ai , 1999]. The second 

improvement expanded the original Ocean Weather 
Station (OWS) P study [Signorini et al . , 2001] to 
include additional state variables (TC0 2 , 0 2 , DOC, and 
Fe) and added a carbonate chemistry sub-model to 
calculate mixed-layer pC0 2 concentrations. The model 
was also upgraded to include air-sea gas exchange 
algorithms featuring five different user-selected 
formulations: Monahan and Spillane [1984], Liss and 
Merlivat [1986], Tans et ai [1990], Wanninkhof [19921 
and Wanninkhof and McGillis [1999]. We will discuss 
the influence of gas exchange formulations on the 
model performance in the forthcoming validation 
section. The fourth, and most recent, version was 
implemented for this study and consists of additional 
state variables (phosphate, dissolved organic carbon, 
dissolved organic nitrogen, and dissolved organic 
phosphorus) specifically implemented in the model for 
the Sargasso Sea application. In addition, we 
implemented a light-dependent Chl:N diagnostic ratio 
(in mg of Chl/mmol N) to account for photoadaptation 
following the parameterization of Doney et ai [1996]: 


Chl:N=Chl:Nr -(CM:Nr <1* 

Chl:N=Ch:N 0 , ^>1* 

where ChliN* 1 ™ =2.0, Chl:N 0 = 1.0, and I*=15 W m 2 . 
These selected values provide the best agreement 
between the BATS chlorophyll record and the model 
prediction. The diagnostic chlorophyll values (Chi = P 
Chl:N) are used to determine the depth attenuation of 

Ipar* 

The model of Gregg and Carder [1990] is used to 
compute the clear sky irradiance in the PAR spectrum 
(350-700 nm). The model is extended to include 
shorter wave lengths (280-700 nm) given the solar 
constants and ozone absorption coefficients for those 
wavelenghts. Cloud cover and meteorological 

parameters required by the spectral surface irradiance 
model were obtained from the NCEP/NCAR reanalysis 
products. The clear sky irradiance is modified to 
account for the observed cloud cover by applying a 
power law correction ( McClain et ai , 1996) tuned to 
yield the observed climatological monthly mean surface 
irradiances ( Dobson and Smith , 1988). Figure 1 shows 
a comparison between model and observed surface 
PAR. The observed values are from 82 cruises taken 
near Bermuda from 1992 to 1997 (D. Siegel personal 
communication). The differences between model and 
observed surface PAR are attributed to small-scale 
cloudiness variability not captured by the coarse 
resolution of the NCEP/NCAR reanalysis products (-2° 
resolution). The water column downwelling irradiance 
is exponentially attenuated by absorption spectra for 
water [Baker and Smith , 1982] and chlorophyll 

[Dupouy et ai , 1997]. 

We also introduced a quadratic formulation [Steele 
and Henderson , 1992; Steele and Henderson , 1995] for 
zooplankton mortality in the model, (g + AZ)Z, which 
provides a more realistic control of the predator 
abundance and phytoplankton concentration at BATS. 

3,1 Governing Equations and Parameterization 

The current model features the following system 
of coupled differential equations to simulate the 
dynamics of phytoplankton nitrogen (P), 
microzooplankton nitrogen (Z), ammonium (NH 4 ), 
nitrate (N0 3 ), phosphate (P0 4 ), dissolved organic 
matter (labile DOC, DON, and DOP), total carbon 
dioxide (TC0 2 ), and oxygen (0 2 ) stocks within the 
upper ocean (surface to 350 meters): 
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Table 1 provides a summary of the model 
variables and input parameter definitions and values 
used for the BATS site. 

The fractions (C/N) Re d, (P/N) Re d, and (0 2 /N) Re d 
follow the Redfield ratios for carbon, nitrogen, 
phosphorus, and oxygen (106C:16N:lP:138O2). The 
symbol N p denotes the net community production. The 
TC0 2 , labile DOC, and oxygen concentration are in 
gravimetric units (pmol kg" 1 ), while all the other 
variables are in volumetric units (mmol m 3 ). The water 
density (p in kg l l ) is calculated using the UNESCO 
International Equation of State (DES 80), as described in 
Fofonoff [1985], which is a function of temperature, 
salinity, and pressure. The Kroenecker delta (S[z = 0] = 
1; 5[z > 0] = 0) is used to denote that the carbon dioxide 
and oxygen fluxes (FC0 2 and F0 2 , respectively) are 
applied at the sea-air interface only. Details of the 
numerical method to solve the coupled differential 
equations are given in McClain et al. [1996]. 
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3.2. Vertical Velocity and Vertical Diffusion 
Formulations 

The vertical velocity variability in the BATS study 
area is affected primarily by vertical displacements of 
the thermocline imparted by the migration of mesoscale 
eddies through the area [McGillicuddy et al . , 1998; 
Siegel et al , 1999]. This eddy pumping mechanism has 
been demonstrated to have a significant impact on the 
nutrient budget at BATS [Siegel et al 1999]. We 
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obtained the eddy pumping vertical velocity time series 
for the BATS location from the divergence of the 
geostrophic velocity field derived from altimeter data. 
Using the thermal wind equations for the geostrophic 
velocity components, combined with the continuity 
equation, we get, 


g 

’ d 

f dHrj ' 

d 

^ dHrj) 

II 

^1 

dx 

l d y ) 

3y 

l dx ). 


where g is the acceleration due to gravity, f is the 
Coriolis parameter, r\ is the SSHA, and H is the 
thermocline depth (-850 m). 

The 10-day high resolution (0.25°) SSHA global 
fields were used to derive the geostrophic currents. The 
SSHA and currents for the mesoscale eddy field around 
Bermuda are shown in Figure 2 for the 10-day average 
periods centered on November 6, 16, 26, and December 
6, 1995. Note that the region near to and north of the 
BATS area is populated by very energetic mesoscale 
eddies. In particular, the BATS site is experiencing the 
effects of a slow moving, anticyclonic/cyclonic eddy 
pair during this time period. This mesoscale eddy 
signature is typical of the region, which is dominated by 
the occurrence of westward propagating features with 
SSHA as large as 25 cm, Eulerian temporal scales of 
roughly a month, lifetimes of several months, spatial 
scales of -200 km, and a propagation of -5 cm s 1 
[Siegel et ai , 1999]. 

Figure 3 shows time series of SSHA and vertical 
velocity (w e dd y ) for the time period of April 1995 
through January 2000, the longest continuous data 
period. The entire data set available at the time covers 
the period of October 22, 1992 - January 14, 2000, but 
with a large gap in the 10-day AVISO composites for 
the period of December 16, 1993 -March 31, 1995, 
during which the SSHA analysis was interrupted in 
favor of the ice and geodesy mission objectives. An 
uninterrupted time series of w ec idy> covering the period 
of 1990-1998, was constructed by filling the missing 
points with data from adjacent years. The vertical 
velocity profile, w(z), is parameterized bilinearly as, 

W(z) = 0.25 w ed dy^’ for z< 200m 

w(z) = w(200) + w eddy - 1 ^ 01 . (13) 

for 200/;; < z< 350m 

This bilinear profile has a break point at 200 m where 
w(z) = 0.25 w c d(iy. At 350 w(z)=w e ddy. The 200 m 
transition depth follows the Siegel et al. [1999] 


empirical parameterization, i.e., below 200 meters, 
isopycnal displacements are related to SLA and above 
200 meters the displacement is modeled as a decreasing 
linear profile that extrapolates to zero at the sea surface. 
The eddy-induced vertical velocity ranges from -3.5 to 
+2.5 m d' 1 (solid line), while the Ekman-induced 
downwelling peaks at -0.1 m d 1 (dotted line). The two 
components of the vertical velocity were added to force 
the model. 

The formulation for the mixed layer eddy 
diffusivity (K v * in m 2 s" 1 ) follows Csanady [1982] with 
an adjusted transitional Reynolds number (25 instead of 
16) for the BATS region. More specifically, 

u * H , , 

K for n*< 0.125 

v 25 

K * , for u*> 0.125 (14) 

V 200 f H mld 


u* = 


where u* is the friction velocity in m s‘ L , x is the wind 
stress in N m' 2 , p w is the water density in kg m \ H m id is 
the mixed layer depth in meters, and f is the Coriolis 
parameter in s' 1 . For numerical stability purposes, the 
value of K v * is bound within the range of 0.0001 and 
0.01 16 in m 2 s' 1 . 

The K v profile is defined as 
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(15) 


where K V b is the maximum of 10% of K v * and 0.0001 
m 2 s' 1 . The parameterizations of w(z) and K v (z) were 
tuned to achieve good agreement between model and 
observed phosphate and nitrate concentrations (below 
detection level near the surface during summer). 


4. Model Forcing and Boundary 
Conditions 

4.1 Physical Forcing 

A continuous time series of physical forcing, 
consisting of daily profiles of ocean variables, daily 
atmospheric variables, and hourly insolation were 
pre-processed and assembled offline to force the 
ecosystem model. Some data sets have time resolutions 
coarser than the time resolution required by the model, 
such that a linear interpolation in time was required to 
fill the gaps. This is the case of the CTD data, which 
are continually obtained at the BATS site on 
approximately monthly intervals. Time series profiles 
of temperature and salinity were interpolated at 1 -meter 
depth intervals and daily time intervals from the 
monthly, 2-meter resolution CTD data sets. Figure 4 
shows the time series of temperature, mixed layer 
depth, eddy diffusivity, vertical velocity, and PAR that 
were used to force the model. The forcing fields were 
produced for the period of 1990-1998 but only 
1992-1998 are shown hereafter because the first two 
years are considered spin-up years and therefore will 
not be analyzed. There are clear interannual changes in 
the forcing fields, most notably in the winter-spring 
mixed layer depths, and accordingly the penetration 
depth of vigorous mixing (K v ), and in the eddy-induced 
upwelling-downwelling events (w e ). 

4.2. Biogeochemical Forcing 

The Neumann boundary condition, 3X/9z=0, is 
applied at both the surface and lower (350 m) model 
domain boundaries for P and Z. Depth-independent 
initial concentrations of P, Z, and NH4 are 0.1, 0.1, and 
0.0 mmol N m'\ respectively. Initial profiles of N0 3 , 
P0 4 , DOC, and DON were obtained from the BATS 
data. The Neumann condition was applied at the 
surface. Since the bottom N0 3 and PO4 concentrations 
are strongly correlated with temperature (r > 0.8), the 
following equations were used to prescribe daily nitrate 
and phosphate concentrations at the bottom boundary 
(350 m): 

NO 3 (350) = 52.315 - 2.676 T(350) (16) 

PO 4 (350) = 3.2084 - 0.168125 T(350) (17) 

For TC0 2 and 0 2 , the observed winter values at 
the surface and at the bottom were used to construct 
initial linear profiles (surface TC0 2 =2050 pmol kg' 1 
and bottom TC0 2 =2090 pmol kg' 1 ; surface O 2 =320 
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pmol kg' 1 and bottom 0 2 =200 pmol kg' 1 ). A fixed 
value equal to the initial condition was applied at the 
lower boundary. The total alkalinity (TA) was 
calculated at each time step as a function of salinity 
following the linear relationship given by Bates et al. 
[1996]: 

TA = -47.155 + 66.576 S, r 2 = 0.91 (18) 

The following formulations for the C0 2 and 0 2 gas 
exchanges were applied in the form of flux boundary 
conditions (FC0 2 and F0 2 in mmol m 2 d 1 ) at the 
sea-air interface: 

FC0 2 =pK 0 aApCO 2 (19) 

F0 2 = pK 0 (0 2 sat -0 2 ) (20) 

where K 0 is the piston velocity, in m d' 1 , which is a 
function of water temperature and wind speed 
(Wanninkhof, 1992; Liss and Merlivat, 1986), a is the 
C0 2 solubility in seawater (in mol kg' 1 atm \ or pmol 
kg' 1 patm 1 ) which is a function of temperature and 
salinity ( Weiss and Price, 1980, Table VI), ApC0 2 (in 
patm) is the difference between air and sea pC0 2 , and 
0 2 sat is the oxygen saturation concentration (in pmol 
kg 1 ) in seawater which is a function of temperature and 
atmospheric pressure (Weiss, 1970). The seawater 
density, p, is used to convert the flux into volumetric 
units. 

The 10-year atmospheric pC0 2 time series used to 
force the model is based on a least-squares fit to the 2.5 
years (mid 1994 to December 1997) of automated 
underway pC0 2 measurements in the vicinity of the 
BATS site (Bates et al, 1998d). The slope and 
intercept of the long-term (secular) trend is a fit to the 
Mauna Loa data with slight adjustments to match the 
interannual pC0 2 trend of the short record at BATS. 
The pC0 2 time series is given by: 


pC O™ = A o + A, t + 

( 2ttt , ^ 
A^sm + Aj 


12 


+ A. sin 
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2 nt 


■ + A c 


( 21 ) 


+A fi sin 
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where t is the time in months (from January 1, 1951), 
Ao = 282.6 patm is the intercept, Ai= 0.125 patm 
month’ 1 is the slope, the amplitudes (A 2 , A 4 , A () ) are 
-7.18, -0.99, and -0.80 patm, respectively. The phases 
(A 3 , A 5 , A 7 ) are 0.86, 0.28, and 0.06 radians, 
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respectively. Figure 5 shows the 1992-1998 time series 
of air pC.0 2 used to force the model, and the 2.5 years 
(1994-1996) of observed daily-averaged air pC0 2 
obtained from underway 30-minute measurements 
(Bates et al, 1998a; Bates et al , 1998b; Bates et al, 
1998c). 

The mixed layer pC0 2 is calculated as a 
function of temperature, salinity, dissolved inorganic 
carbon, and alkalinity following the carbonate 
chemistry approach outlined in Signorini et al. [2001]. 
The TC0 2 concentration depends on the balance 
between air-sea C0 2 flux, the biological uptake of 
carbon, and the vertical influx of TC0 2 through the 
bottom boundary. The biological uptake of carbon is 
performed via the net community production, N p , which 
consists of the sum of all sources and sinks in the 
ammonium and nitrate pools. From (4) and (5) we can 
write 


N p =GP-(a p m + r p )P- 

[ a z (g + hZ) + r z +c pel y /]z- 


k rc DOC 


( 

V C ;Red 


( 22 ) 


Equation (22) is applied to the TCO 2 and O 2 equations 
((7) and (8)), after multiplication by the corresponding 
Redfield ratios, as a biological sink and source of DIC 
and oxygen, respectively. 


5. Model Validation 


5.1 Gas Exchange Formulation 


model and data ApC0 2 and TCO 2 ranges with error in 
percent, and the bias and root-mean-square (rms) error 
for the TC0 2 model prediction. The results of Table 2 
show that the Wanninkhof [1992] algorithm provides 
the smallest errors for the prediction of ApC0 2 range 
and the least TC0 2 bias (+2.5% and -0.01 pmol kg' 1 , 
respectively). The Monahan and Spillane [1984] 
algorithm does not match the observed pC0 2 and TC02 
ranges as well as the other four. In view of these results, 
we adopted the Wanninkhof [1992] algorithm for the 
air-sea gas transfer in the model. Figure 7 shows time 
series of predicted and observed TC0 2 and ocean pC0 2 . 
There is very good agreement with the observed values. 
Monthly time series of ApC0 2 and FC0 2 for 1995-1996 
are shown in Figure 8. The black dots connected with 
solid lines are based on observations. The ApC02 values 
are based on underway air and ocean pCC >2 
measurements (Bates et al., 1998a; Bates et al, 1998b; 
Bates et al., 1998c). The FC0 2 is calculated from the 
observed ApC0 2 and gas transfer coefficient from the 
model simulation. The dashed lines are the ApC0 2 and 
FCO 2 model predictions for the same time period. 
Phase and amplitude of model predictions agree well 
with the observed values. To compare model and 
observed mean values for 1995-1996, one must first 
eliminate January and December 1995, and January 
1996, from the mean, as there are no observed ApCC >2 
values for these 3 months. After doing that, the ApC0 2 
and FC0 2 model means for 1995-1996 are -3.8 patm 
and -0.20 mol C m' 2 yr' 1 , respectively. These values are 
within 16% of the observed ApC0 2 , and within 20% of 
the FC0 2 derived from the observed ApC0 2 and the gas 
transfer coefficient used in the model (-4.5 patm and 
-0.25 mol C m‘ 2 yr' 1 , respectively). 


As stated in Section 3, sensitivity tests were 
conducted using five different gas exchange 
formulations. Figure 6 provides a comparison of 5 
published air-sea gas transfer formulations. The piston 
velocity was converted to m yr' 1 to conform with the 
original units (mol m' 2 yr' 1 patm' 1 ) of the Tans et al. 
[1990] gas transfer coefficient, which is independent of 
solubility and Schmidt number. The comparison is 
made at constant salinity (36.0) and temperature (20°C). 
The behavior of the algorithms is quite different, which 
is exacerbated at wind speeds greater than 10 m s' 1 . 
However, within the range of 4 to 11 ms' 1 , the Tans et 
al. [1990] and Wanninkhof [1992] algorithms are very 
similar. 

Table 2 summarizes the results of the five 
sensitivity tests using the gas exchange formulations 
shown in Figure 6. The table shows the mean value and 
range of the sea-air flux (FCO 2 ), a comparison between 


5.2. Validation of Biogeochemical Parameters 

Seasonal climatological profiles of predicted and 
observed biogeochemical parameters are shown in 
Figures 9 through 17. The observed data points (black 
dots) covering the period of 1988-1998 are sorted by 
season and averaged into seasonal profiles (dotted 
lines). The predicted seasonal climatological profiles 
are shown as solid lines. Inspection of Figures 9 
through 17, reveals that all biogeochemical variables 
are in very good agreement with model predictions. 

Nutrients at the BATS site are below colorimetric 
detection limits [Bates et al, 1996] most of the year 
(< 0.05 pmol kg' 1 ). This nutrient depletion is also 
manifested in the model predictions as shown in the 
seasonal profiles of nitrate and phosphate (Figures 9 
and 10). The predicted chlorophyll profiles have 
surface concentrations and subsurface maximum values 
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in winter and summer that agree well with the 
equivalent features exhibited by the data. These 
features are less comparable with the data in summer 
and fall where the model tends to provide higher surface 
chlorophyll values and shallower chlorophyll maxima. 
The predicted seasonal profiles of primary production, 
TC0 2 , 0 2 , and DON are all within the bounds of the 
data. Table 3 summarizes the comparison between 
predicted and observed climatological averages of 
surface PAR and biogeochemical parameters separated 
into seasons and yearly averages. The predicted surface 
PAR is within 7% of the observed value. The model 
overpredicts the surface chlorophyll by 30%, mostly 
due to the summer-fall discrepancies. The prediction 
accuracy of the nitrate and phosphate concentrations 
(0-50 m) are hard to evaluate because the observed 
values are biased due to inadequate sensitivity. The 
annual model nitrate (0-50m and 1992-1998 averages) 
is almost 3 times the observed value, while the 
phosphate is within 13% of the observed value. The 
predicted primary production (0- 1 50 m), surface TC0 2 , 
surface 0 2 , surface DON, and surface DOC are within 
7%, 0.2%, 0.8%, 8%, and 1.4% of the observed values, 
respectively. 


6. Discussion 

6.1 Seasonal and Interannual Variability 

The seasonal variability and interannual variability 
of the biogeochemical parameters respond to the 
variability of the physical forcing. Figure 18 highlights 
the biogeochemical-physical interactions on a seasonal 
and interannual basis. This figure shows a composite 
time series of PAR at 50 m, mixed-layer depth (MLD), 
SST, surface TC0 2 , surface TA, ApC0 2 , air-sea C0 2 
flux (FC0 2 ), net community production (N p , negative 
sign denotes a drawdown in C0 2 ), new production, 
primary production, and N0 3 and P0 4 concentrations 
(0-50 m). Except for MLD, SST, and TA, all of the 
above parameters were predicted by the model. As 
mentioned in Section 4, TA is calculated following 
Bates et al. [1996], and the MLD and SST are derived 
from Hydrostation S data. 

The long-term patterns of MLD and mixed-layer 
temperature at BATS show significant interannual 
fluctuations and coherent variations on longer time 
scales. Some of the variations in MLD also seem 
correlated with the El Nino-Southern Oscillation 
[Michaels and Knap, 1996], The years of deep mixed 
layers in the first half of the record (1955-1980) at 
Hydrostation S are often coincident with an El Nino 
event and the years of extremely shallow mixed layers 


coincide with the cold events (La Nina). After 1980, 
this pattern seems less clear but coherent with the 
change in patterns of El Ninos [Michaels and Knap, 
1996]. The depth and temperature of the winter MLD 
play a major role in determining interannual variations 
in total primary and new production [Menzel and 
Ryther, 1961; Michaels et al, 1994a, Michaels et al., 
1994b]. Figure 18 shows that, in years of deep and cold 
mixed layers, the mixed-layer nitrate and phosphate 
concentrations are significantly elevated above the 
concentrations in years with shallow and warm mixed 
layers. This interannual scenario agrees well with the 
data [Michaels et al., 1994a]. However, since all of this 
near-surface nitrate and phosphorus are quickly 
consumed by phytoplankton, their uptake represents a 
minimum amount of new production. The net primary 
production is significantly larger, indicating that most 
of the required nutrients are recycled by the community. 
Thus the model is configured to recycle 80% of the 
required nutrients from organic matter released via 
zooplankton and phytoplankton mortality and other 
organic waste (zooplankton fecal pellets and excretion 
and phytoplankton extracellular release). The other 
20% of organic matter is equally divided into DOM and 
POM exported from the euphotic zone. These fractions 
seem to be adequate judging from the good agreement 
between model and observed nutrient concentrations 
and uptake rates. 

A strong seasonality is evident in the in the sea-air 
ApC0 2 , with negative values in winter when SST is 
minimum, and positive in summer when the MLD is 
shallow and SST is warmest. The seasonal cycle of 
FC0 2 has an amplitude of -5 mmol C m 2 d 1 with peak 
outgassing in summer and peak ingassing in winter. 
The net community production has a similar seasonal 
amplitude but out of phase with FC0 2 . Production 
peaks in spring and reaches a minimum in the fall. 

The seasonal variability of surface layer TA is 
quite large. For the period of model simulation 
(1992-1998), the TA ranged between a low of 2364 
pmol kg' 1 and a high of 2404 pmol kg' 1 . A significant 
portion of the TA variability is due to salinity changes 
via the effects of evaporation and precipitation [Bates et 
al, 1996]. The minimum TA (2364 pmol kg' 1 ) 
occurred on August 1998, during the peak of La Nina. 
This minimum TA was accompanied by the warmest 
SST (28.4°C), lowest surface salinity (36.23), highest 
pC0 2 (420 patm), and lowest TC0 2 (2018 pmol kg' 1 ) of 
the period. This low TC0 2 concentration resulted from 
the combination of warmest SST and lowest TA, 
producing the minimum dissolved C0 2 concentration of 
the period. 

Time series of vertical profiles of chlorophyll 
concentration, primary production, phosphate, nitrate, 



zooplankton, ammonium, oxygen and total carbon 
dioxide are shown in Figures 19 and 20. The vertical 
stratification of these parameters portrays a classical 
scenario of an oligotrophic environment. Nutrients are 
depleted within the euphotic zone most of the time and 
there is a distinct sub-surface chlorophyll maximum, 
which peaks in summer-fall. There is also a seasonal 
pattern in the ammonium vertical stratification with 
sub-surface (125 m) peaks during March when the 
mixed layer is deepest and primary production reaches 
its peak. The euphotic zone oxygen concentration 
peaks in February-March when solubility is highest 
(lowest temperatures) and primary production is also 
highest. The seasonal pattern of TC0 2 is characterized 
by minimum concentration within the upper 50 m 
during August when the mixed layer is shallowest and 
near-surface temperatures reach their peak. The TC0 2 
concentration increases monotonically towards the 
bottom (350 m) where it is set, via the boundary 
condition, to a concentration of -2090 pmol kg' 1 . 
Below the euphotic zone, the phosphate and nitrate 
concentrations exhibit an interannual pattern coherent 
with the physical forcing (eddy upwelling and mixing). 
The phosphate and nitrate below -100 meters shows 
interannual changes that depend upon the degree of 
winter mixing (e.g., mixed layer depth), nutricline 
vertical displacement, and the accompanying 
upwelling/downwelling induced by the mesoscale 
eddies. 

6.2. Nitrogen and Phosphorus Budgets 

The climatological (1992-1998) nitrogen and 
phosphorus fluxes predicted by the model are illustrated 
in the block diagrams in Figures 21 and 22. This 
diagram shows the pathways within the ecosystem and 
the exchanges through the bottom boundary (350 m). 
The mean values and standard deviations (in mmol m 2 
yr' 1 ) are shown for each of the nitrogen and phosphorus 
pathways. The nutrient supply at the bottom boundary 
is divided into a diffusive/convective component and an 
advection component. The effects of convective and 
diffusive mixing are combined because the vertical 
eddy diffusivity in the model effectively takes into 
account both mixing processes. These components 
have the highest interannual variability since their 
standard deviations are close to or larger than the mean. 
The upward nutrient flux through the bottom boundary 
is balanced by sedimentation, 334±76 mmol m' 2 yr* 1 for 
nitrate and 21±8 mmol m' 2 yr 1 for phosphate. 
Approximately 30% of the required nutrients are 
transported upward via advection, the other 70% are 
provided by convective and diffusive processes. 
Convection of nutrients reaches a maximum during 
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winter-spring when the mixed layer is deepest. 
Separation of the eddy-upwelling component in the 
model is difficult because vertical excursions of the 
nutricline from the passage of mesoscale eddies affect 
both advective and diffusive components of the nutrient 
flux. Therefore, the effect of eddy upwelling may be 
larger than 30% (advective component) and approach 
the -50% reported by Siegel et al [1999]. The new 
production predicted by the model, shown in Figure 21, 
nitrogen uptake of 0.47±0.14 mol N m' 2 yr' 1 , is within 
6% of the 0.50±0.15 mol N m' 2 yr* 1 mean estimate for 
the Sargasso Sea reported by Siegel et al [1999]. The 
predicted and estimated values are statistically identical 
judging from the magnitude of their individual standard 
deviations. 

6.3. Carbon Budget 

Figure 23 shows the climatological carbon balance 
for the BATS region based on the 7-year (1992-1998) 
model simulation. The two boxes in Figure 23 
represent the dissolved inorganic carbon (DIC) in 
surface water (solubility pump), and the biological 
pump, which takes up C0 2 from the DIC pool via 
photosynthesis and exports POC and labile DOC as a 
byproduct. The labile DOC is quickly consumed by 
bacteria and remineralized into DIC. The biological 
DIC consumption, i.e., the net community production 
(2.3 mol C m 2 yr' 1 ), is balanced by the air-sea C0 2 flux 
(0.3 mol C m* 2 yr* 1 ), the remineralization of DOC (0.9 
mol C m' 2 yr' 1 ), and the DIC flux (1.1 mol C m' 2 yr' 1 ) at 
the bottom boundary (350 in). In terms of overall 
carbon balance, the sedimentation rate (POC + fecal 
pellets), 1.4 mol C m' 2 yr* ! , is balanced by the air-sea 
C0 2 flux and the DIC advection through the bottom 
boundary. Michaels et al. [1994] calculate a POC flux 
of 0.4 to 1.9 mol C m* 2 for the period April -December, 
based on deficits of 234 Th. This POC flux estimate is 
for 150 meters and it is based on 1988-1993 BATS data. 
The model POC flux (1992-1998) for the same depth 
and time-averaging constraints (April-December) 
ranges from 0.6 to 0.8 mol C m* 2 , which is within the 
range of the Michaels et al [1994] estimate. 

7. Summary and Conclusions 

An analysis of biogeochemical- physical 
interactions and carbon cycle was conducted for the 
BATS region using a combination of field observations 
and a one-dimensional ecosystem/carbon-cycle model. 
The required forcing and boundary conditions were 
derived from the BATS comprehensive hydrographic 
and biogeochemical data sets. The model performance 
was evaluated by comparing 10 predicted parameters 
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(ApC0 2 , TC0 2 , 0 2 , PP, PAR, Chl-a, N0 3 , P0 4 , DON, 
DOC) with available time series of equivalent 
parameters observed at BATS. Most model results fall 
within 8% of the observed values. 

Primary production peaks at 657 mg C m' d' in 
spring and is reduced to 375 mg C m' 2 d" 1 in the fall. 
The surface nitrate and phosphate concentrations are 
extremely low during summer-fall (0. 1-0.2 and 
0.002-0.003 mmol m'\ respectively). Regenerated 
production is set to a very high level in the model 
(-80% of the nutrients are recycled) in order to achieve 
a good match with observations, an indication of the 
oligotrophic nature of the BATS site. The model also 
shows that a first-order nitrogen balance can be 
achieved without resorting to nitrification processes 
from other species, such as cyanobacteria nitrogen 
fixation. 

The model results reveal significant interannual 
variability of mixed-layer TC0 2 . There was an upward 
trend of -1.2 pmol kg' 1 yr' 1 in the surface TC0 2 during 
1992-1996, which is similar to the -1.7 pmol kg' 1 yr' 1 
observed by Bates et al. [1996] during 1988-1993. 
They attributed this trend to the uptake of atmospheric 
C0 2 through gas exchange, or natural variability of the 
subtropical gyre. However, our model results show that 
this trend is reversed during 1997-1998 to -2.7 pmol 


kg' 1 yr 1 . This relatively large TC0 2 change was 
accompanied by a significant drop (9 pmol kg ') in TA 
from 2392 pmol kg' 1 in 1997 to 2383 pmol kg' 1 in 
1998. The minimum TA (2364 pmol kg 1 ) occurred on 
August 1998, during the peak of La Nina. The drop in 
TA was large enough to reverse the TC0 2 trend 
imposed by atmospheric C0 2 uptake. 

The sea-air C0 2 flux was evaluated using 5 
different gas transfer coefficients. The C0 2 flux 
obtained with these 5 methods ranges from -0.12 to 
-0.50 mol C m 2 yr' 1 , indicating that the BATS region is 
a weak sink for atmospheric C0 2 . The mean C0 2 flux 
for 1992-1998, -0.3 mol C m' 2 yr 1 , was calculated using 
the Wanninkhof [ 1992] gas transfer coefficient because 
it provides the best match of ApC0 2 with observations 
and the least TCO z bias (model-observation annual 
means). The atmospheric uptake of C0 2 is balanced by 
an upward bottom (350 m) flux of DIC (+1.1 mol C m 
yr 1 ) and sedimentation (-1.4 mol C m' 2 yr 1 ). These 
results show that it is possible to achieve a first-order 
carbon cycle balance without considering the effects of 
horizontal advection. 
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Table 1 . Model variables and input parameters. For variables (P, Z, NH 4 , NO 3 , DON, DOC, DOP, TC0 2 , and 
O 2 ), the column of values contains initial mixed-layer concentrations values (for DON, DOC, DOP, TCO 2 , and 
O 2 , bottom concentration values are also provided). The remaining values are the values assigned to the various 
model parameters. 


Symbol Value Definition 


p 

0 . 1 / 0.1 

Phytoplankton concentration (mmol N m -3 ) 

z 

0 . 1 / 0.1 

Zooplankton concentration (mmol N m“ 3 ) 

NHL, 

o.o/o.o 

Ammonium concentration (mmol m -3 ) 

N0 3 

0.2/4.0 

Nitrate concentration (mmol m -3 ) 

P0 4 

0 . 01 / 0.20 

Phosphate concentration (mmol m -3 ) 

DON 

1 . 0 / 0.0 

Labile dissolved organic nitrogen (mmol m~ 3 ) 

DOC 

17.0/0.0 

Labile dissolved organic carbon (//mol kg^ 1 ) 

DOP 

0.039/0.0 

Labile dissolved organic phosphorous (mmol m -3 ) 

TC0 2 

2050/2090 

Total dissolved inorganic carbon (//mol kg -1 ) 
Oxygen concentration (//mol kg -1 ) 

0 2 

320/200 

TTl 


Regenerated production fraction, Dimensionless 

7T 2 


New production production fraction, Dimensionless 

KnO 3 

0.024 

Half saturation for nitrate uptake (mmol m“ 3 ) 

KnH 4 

0.024 

Half saturation for ammonium uptake (mmol m -3 ) 

m 

0.05 

Phytoplankton mortality rate (d -1 ) 

Go 

0.5899 

Phytoplankton growth rate at CPC (d -1 ) 

k 

' '■gp 

0.0633 

Temperature sensitivity of algal growth (°C -1 ) 

pk 

5 

Ammonium inhibition of nitrate uptake, Dimensionless 

c 

0 

Maximum phytoplankton sinking rate (m d -1 ) 

9 

0.12 

Linear coefficient of zooplankton mortality (d^ 1 ) 

h 

0.30 

Quadratic coefficient of zooplankton mortality (mmol ' 1 m 3 d _1 ) 

1 


Herbivore grazing term, I = i^AP[l - exp(—AP)\ } (d” 1 ) 
Maximum grazing rate (d _1 ) 

R 

7.0 

A 

2.0 

Herbivore grazing constant (m 3 mmol N -1 ) 

7 

0.3 

Unassimilated zooplankton ingestion ratio, Dimensionless 

T ZO 

0.019 

Respiration rate for zooplankton at 0°C (d' 1 ) 

k rz 

0.08 

Temperature sensitivity of zooplankton respiration ( 0 C _1 ) 

dp 

0.8 

Fraction of dead phytoplankton converted to ammonium, Dimensionless 

dz 

0.8 

Fraction of dead zooplankton converted to ammonium, Dimensionless 

a' p 

0.1 

Fraction of dead phytoplankton converted to DOC, DON, and DOP, Dimensionless 

a'z 

0.1 

Fraction of dead zooplankton converted to DOC, DON, and DOP, Dimensionless 

k rc 

0.00005 

Remineralization rate of DOC (d _1 ) 

1 » 

A >rn 

0.0001 

Remineralization rate of DON (d _1 ) 

k 

n, r p 

0.0001 

Remineralization rate of DOP (d _1 ) 

An 

Gi max 

2.0 

Maximum rate of nitrification (mmol N d _1 ) 


0.0095 

Minimum inhibition dosage for nitrification (W m“ 2 ) 

k d 

0.036 

Half saturation dosage for nitrification photoinhibition (W m~ 2 ) 

dpel 

0.8 

Fecal pellet remineralization fraction, Dimensionless 
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Table 2. Mean sea-air C0 2 flux (FC0 2 ), mean range of sea-air C0 2 flux, and error analysis for model-predicted 
ocean ApC0 2 and TC0 2 using five gas transfer coefficients (K) and BATS site observations as a reference. Data 
availability limits the ApC0 2 model-data comparison to 1995-1996 and the TC0 2 model-data comparison to 
1992-1997. The TC0 2 bias is defined as the difference between the model and observed means for 1992-1998. 
The RMS error is also calculated for 1992-1998. 



FC0 2 , 

mol m 2 yr 1 

ApCO ; 

2 Range, fi atm 

tco 2 

Range, /xmol kg 1 

tco 2 , 

^mol kg 1 

Mean 

Range 

Model 

Data 

Error 

Model 

Data 

Error 

Bias 

RMS 

K \ 

-0.50 

3.93 

91.7 

94.2 

- 2 . 7 % 

37.8 

46.2 

-18% 

+0.62 

±8.0 

k 2 

-0.44 

3.40 

96.6 

94.2 

+ 2 . 5 % 

33.5 

46.2 

-28% 

-0.01 

±8.0 

k 3 

-0.32 

2.12 

104.8 

94.2 

+ 11 . 2 % 

26.4 

46.2 

-43% 

-1.02 

±8.9 

k 4 

-0.40 

3.17 

103.8 

94.2 

4 - 10 . 1 % 

28.0 

46.2 

-39% 

+0.85 

±8.8 

K 5 

-0.12 

0.79 

107.7 

94.2 

4-14.3% 

21.2 

46.2 

-54% 

-5.48 

±11.6 


The gas exchange formulations used in the analysis are: K\=^Tans et al [1990]; 1 ^ 2 = Wanninkhof [1992]; Kz—Liss 
and Merlivat [1986]; K 4 = Wanninkhof and McGillis [1999]; and K^ = Monahan and Spillane [1984]. 
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Table 3. Summary of seasonal and climatological comparisons between predicted and observed biogeo chemical 
parameters. The climatological values from the data are based on 1988-1998 averages, while those obtained from 
the model are based on 1992-1998 averages. 



Winter 

Spring 

Summer 

Fall 

Year 

Model 

Data 

Model 

Data 

Model 

Data 

Model 

Data 

Model 

Data 

PAR 

46.0 

49.3 

83.5 

89.8 

96.3 

99.2 

61.5 

68.9 

71.8 

76.8 

Chl-a 

0.14 

0.16 

0.14 

0.11 

0.11 

0.04 

0.12 

0.07 

0.13 

0.10 

N0 3 

0.43 

0.16 

0.64 

0.11 

0.16 

0.06 

0.09 

0.16 

0.33 

0.12 

P0 4 

0.012 

0.009 

0.017 

0.003 

0.002 

0.010 

0.003 

0.008 

0.009 

0.008 

PP 

456.3 

432.4 

656.5 

593.5 

440.4 

381.8 

375.0 

387.3 

482.1 

448.8 

DON 

0.94 

0.94 

0.88 

1.05 

0.97 

1.10 

1.00 

1.03 

0.95 

1.03 

DOP 

0.06 


0.05 


0.06 


0.06 


0.06 


DOC 

14.0 

12.9 

13.1 

13.1 

14.4 

15.3 

14.4 

15.5 

14.0 

14.2 

TC0 2 

2049.4 

2046.2 

2058.1 

2055.3 

2036.9 

2035.5 

2030.4 

2024.1 

2043.7 

2040.3 

O2 

223.2 

218.0 

223.7 

223.8 

207.2 

206.0 

208.3 

206.9 

215.6 

213.7 


Nitrate, phosphate, labile DON, DOP (mmol/m 3 ), and DOC (/nnol kg -1 ) were averaged from surface to 50 m. The 
Chl-a (mg/m 3 ), TC02, O 2 (/i mol kg” 1 ) concentrations are surface values. The observed surface irradiance (PAR[0 - ]), in 
W m~ 2 , is obtained from 1992-1997 averages. The primary production (PP in mg C m -2 d~*) is integrated from surface 
to 150 m for the period of 1992-1997. 



PAR(O-), W/m 2 
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1992 1993 1994 1995 1996 1997 1998 


Figure 1: Comparison between model and observed daily averaged surface PAR (just below the surface at 
m) for 1992-1997. The solid line is the model time series while the black circles connected by the dotted 
line are the observed values. 
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Figure 2: Sea surface height anomaly (cm) and geostrophically derived currents near the BATS site for 
four 10-day composites in November-December 1995. Note the two mesoscale eddies, one cyclonic and one 
anticyclonic, fringing the BATS site (white circle). 
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Figure 3: Time series of 
vertical velocities (bottorr 
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Figure 4: Time series of temperature (T), mixed layer depth (white line), vertical eddy diffusivity (ivy) 
vertical velocity (tu), and PAR for the period of 1992-1998. 




19 


Model (solid) and Observed (black circles) Air p0O 2 



Figure 5: Time series of air pC0 2 for 1992-1998. The thick line denotes the time series used to force the 
model. The black circles are daily-averaged air pC0 2 values obtained from underway 30-min measurements 
(Bates et al , 1998a, 1998b, 1998c). 
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* K = a 10.33 W 3 (Sc/Sc20)*° 5 , Wanninkhof and McGlllis [1999] 

K = a 27.16 W 2 (Sc/Sc20)-° * Wanninkhof [1992] 

K = 0.016 (W - 3), Tana et at. [1990] 

K = a (839.5 + 0.456 W 3 ) (Sc/Sc20)*°- 5 , Monahan and Splllane [1984] 
Lisa and Meriivat [1986] 

K = a 1 4.892 W (Sc/Sc20)^ 5 , for W <; 3.6 
K = a (249.66 W - 845.34) (ScySc20)-°* 5 , for 3.6 < W ^ 13 
K = a (516.84 W - 4318.68) (Sc/Sc20)-°- 5 , for W > 13 



T 

10 

Wind Speed (m/s) 


Figure 6: Gas transfer coefficient as a function of wind speed. Curves are shown for five different formulations. 
All constants have been converted to reflect units of mol C0 2 m -2 yr -1 /iatm~\ which are the units of the 
Tans et al. [1990] formulation. The assumptions are: SST=20°C and SSS=36, so that the C0 2 solubility 
a=3.223xl0~° mol m -3 /iatm -1 and the Schmidt number ratio (S/Sc20)~° 5 =1.0. 
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Model (Solid) and Observed (dotted) Surface TC0 2 
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Figure 7: Time series of TC0 2 and ocean pC0 2 for 1992-1998. The thick lines denote model predictions. 
The observed TC0 2 (1992-1997) is shown by the dotted line. The black circles are daily-averaged ocean 
pC0 2 values obtained from underway 30-min measurements ( Bates et al., 1998a. 1998b, 1998c). 
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Observed = black circles Model = dashed line 



Figure 8: Monthly time series of ApC0 2 (top panel) and sea-air C0 2 flux (FC0 2 , bottom panel) for 1995- 
1996. The black dots connected with solid lines are based on observations. The ApC0 2 values are based on 
underway air and ocean pC0 2 measurements ( Bates et al , 1998). The FC0 2 is calculated from the observed 
ApC0 2 and gas transfer coefficient from the model simulation. The dashed lines are the model predictions 
for the same time period. 
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Winter DON (|iM) Spring DON (jiM) 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 15: Same as in Figure 9, except for winter, spring, summer, and fall dissolved organic nitrogen (DON). 
Since there are no direct observations of DON, the DON ” observed” seasonal profiles were obtained from 
the difference between the TON and the PON profiles. Thus no data points are shown, just the averaged 
profile (dotted line). Labile DON was obtained by subtracting the total DON observed mean value at 350 
m. 








Winter OOP (pM) Spring DOP (pM) 
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Summer DOP (pM) Fall DOP (pM) 
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Figure 16: Same as in Figure 15, except for winter, spring, summer, and fall dissolved organic phosphate 
(DOP). There are no observed values of DOP. 
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Figure 18: Time series of daily predicted PAR at 50 m, observed mixed layer depth (MLD, green), sea 
surface temperature (SST, red), surface TCO 2 (blue), surface TA (green), ApC0 2) sea-air C0 2 flux (FC 02 , 
red line), net community production ( N P) blue line), new production (NewP, green line), net primary 
production ( NetP , blue line) for the whole water column (blue line), and phosphate (PO 4 ) and nitrate 
(NO 3 ) at 50 m. 
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Figure 19: Monthly time series profiles of predicted Chl-a, primary production (PP), phosphate (PO 4 ), and 
nitrate (N0 3 ) for 1992-1998. 
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Figure 20: Same as in Figure 19, except for predicted zooplankton (Z), ammonium (NH 4 ), oxygen (0 2 ), 
and total carbon dioxide (TCOo). 
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Figure 22: Climatological phosphate flux (mmol P m'^ yr"*) for the BATS study 
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